Transport theory yields renormalization group equations 



Jochen Rau 

Max- Planck- Institut fiir Physik komplexer Systeme, Bayreuther Strafie 40 Haus 16, 01187 Dresden, Germany 
(revised version, January 21, 1997; to appear in Phys. Rev. E) 



OO 



I 

o 
o 



(N 

> 

On 



X3 
O 

o 



X 
J3 



We show that dissipative transport and renormalization 
can be described in a single theoretical framework. The ap- 
propriate mathematical tool is the Nakajima-Zwanzig projec- 
tion technique. We illustrate our result in the case of inter- 
acting quantum gases, where we use the Nakajima-Zwanzig 
approach to investigate the renormalization group flow of the 
effective two-body interaction. 
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I. INTRODUCTION 

The basic theme of statistical mechanics -how to ob- 
tain a system's macroscopic properties from the laws of 
its underlying microscopic dynamics- appears in many 
variations. Two out of many examples are the prob- 
lem of determining critical exponents at second-order 
phase transitions, or the problem of deriving macroscopic 
transport equations. The former is usually tackled with 
the help of Wilson's renormalization group,llrQ a math- 
ematical tool that allows one to iteratively eliminate 
short-wavelength modes and thus to arrive at effective 
( "renormalized" ) theories which describe the dynamics 
on successively larger length scales. The latter has been 
tackled in various ways, among [them the so-called pro- 
jection technique by Nakajima,El Zwanzigjj Morj3 and 
Robertson.! Eliminating unmonitored, rapidly oscillat- 
ing degrees of freedom from the equation of motion by 
means of suitable projections in the space of observables, 
the projection technique yields closed (but generally no 
longer Markovian) "transport equations" for the selected 
macroscopic degrees of freedom. 

While the two methods -Wilson's renormalization 
group and Zwanzig's projection technique- may appear 
quite different in their mathematical formulation, they 
are very similar in spirit. In both cases one strives to 
focus on selected features of the dynamics (its infrared 
limit, or the evolution of only few macroscopic observ- 
ables) deemed interesting, and to this end devises a sys- 
tematic procedure for eliminating all other, "irrelevant" 
degrees of freedom (a procedure commonly referred to as 
"coarse graining"). Discarding thus unnecessary "bag- 
gage" from the problem at hand, one succeeds in describ- 
ing the interesting features of the dynamics without ever 
having to solve the complete, and far too complicated, 
microscopic theory. This similarity of the basic approach 



suggests that renormalization and the transition from mi- 
crodynamics to macroscopic transport are in fact closely 
related procedures, and that it should be possible to cast 
them into a common theoretical framework. 

Building a bridge between renormalization and trans- 
port theory would be not only satisfying conceptually, 
but would also help tackle a variety of practical problems. 
Often, the macroscopic evolution of a complex quantum 
system exhibits both dissipation and modified, "renor- 
malized" dynamical parameters such as effective masses 
or effective interactions. Let us consider, for example, 
liquid ■^He or nuclear matter away from equilibrium. In 
order to formulate a macroscopic transport theory for 
such an interacting fermion system one must perform two 
consecutive coarse graining procedures: first eliminating 
short- wavelength modes to arrive at an effective ( "renor- 
malized" ) theory for quasiparticle excitations close to the 
Fermi surface, which typically feature effective masses 
and screened interactionspllJ and then discarding their 
statistical correlations to obtain an Uehling-Uhlenbeck- 
type transport equation for the single-quasiparticle dis- 
tribution. The description of the macroscopic dynamics, 
therefore, requires an appropriate combination of renor- 
malization and statistical coarse graining. 

Clearly, the two coarse grainings do not commute; the 
latter (statistical coarse graining) is contingent upon the 
former (renormalization). For instance, the renormal- 
ization group flow yields screeningtil and hence renders 
the interaction range finite, thus generating that separa- 
tion of scales which is indispensable for thc-subscquent 
derivation of a Markovian transport theory£3 But what 
happens if scales converge rather than separate upon 
renormalization? How are then renormalization and sta- 
tistical coarse graining best combined? More generally, 
what is the connection between effective dynamics and 
dissipation? Does their interplay lead to interesting new 
phenomena? To what extent can renormalization group 
techniques be applied to study nonequilibrium, dissipa- 
tive processes? How do transport coefficients change 
under renormalization group transformations? Is it al- 
ways true that transport coefficients are renormalized 
by simply trading bare masses and couplings for their 
renormalized counterparts, while keeping the £arm of 
the functional dependence on these parameters?E3 And, 
somewhat speculative, are there "universality classes" of 
transport theories? These and other issues might be best 
approached in a unified mathematical framework that 
encompasses both renormalization and dissipative trans- 
port as special cases. 

There has already been some progress towards such 
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a unified picture. The success of Andersonis, "poor 
man's scaling" approach to the Kondo problem,0 Seke's 
projeclion-method treatment of the nonrelativistic Lamb 
shift ,E3 the calculation of the 1-loop renormalization of 
0* theory by means of Bloch-Fcshbach techniques or a 
recent renormalization group study of interacting fermion 
systems (BCS instability, screening) within a purely alge- 
braic framework ,EZl suggest that one can formulate Wil- 
son's renormalization in terms of projections in Hilbert 
space, completely analogous to the projections in the 
space of observables which, in the Nakajima-Zwanzig ap- 
proach, would lead to macroscopic transport equations. 

In the present paper I wish to make this analogy even 
more explicit. I will show that one can actually obtain 
renormalization group equations within the Nakajima- 
Zwanzig projection approach, and that hence renormal- 
ization can be embedded into the general mathematical 
framework of transport theory. After a brief introduc- 
tion to the projection technique (Sec. f|) and a discus- 
sion of various approximations (Sec. Ol) I shall isolate 
the dissipative and non-dissipative parts of the macro- 
scopic dynamics, and show that the latter is governed 
by an effective, "renormalized" Hamiltonian (Sec. IV). 
For illustration, these general ideas are then applied to 
studying the low-energy dynamics of interacting quan- 
tum (Bose and Fermi) gases, in particular the renormal- 
ization group flow of their effective two-body interaction 
(Sees. and VI ) . Fi nally, I shall conclude with a brief 
summary in Sec. VII. 



II. PROJECTION TECHNIQUE 

In this section I give a very brief intcoduction to 
the Nakajima-Zwanzig projection technia«e.Eril More de- 
tails canu-he, found in several textbookstS and in recent 
reviews. E3o 

When studying the dynamics of a macroscopic quan- 
tum system away from equilibrium, one typically moni- 
tors the evolution of the expectation values 



ga{t) tv[pit)Ga] 



(1) 



of only a very small set of selected ( "relevant" ) observ- 
ables {Ga}- These evolve according to 

gait) ^ iipit)\CGa) , (2) 

with p(t) being the statistical operator, C the Liouvillian 

C:=n-^[H,*] (3) 

associated with the Hamiltonian and the inner prod- 
uct ('I') defined as 



{A\B) := tT[A^B] 



(4) 



The equation of motion in the form (g) does not con- 
stitute a closed system of differential equations for the 



selected expectation values {ga{t)}; its right-hand side 
will generally depend not just on the selected, but also 
on all the other ("irrelevant") degrees of freedom. With 
the help of the projection technique to be sketched be- 
low, these irrelevant degrees of freedom can be systemati- 
cally eliminated from the equation of motion, in exchange 
for non-Markovian and (possibly) nonlinear features of 
the resulting closed "transport equation" for the {ga{t)}- 
Mapping thus the influence of irrelevant degrees of free- 
dom onto -among other features- a non-local behavior in 
time opens the way to the exploitation of well-separated 
time scales, and hence serves as a good starting point 
for powerful approximations such as the Markovian and 
quasistationary limits. Indeed, in this fashion one can 
derive many of the well-known equations of nonequilib- 
rium statistical mechanics, for example rate equations, 
the quantum Boltzmann, Master, Langevin-Mori or even 
time-dependent Hartree-Fock equations. And as I will 
show later in this paper, the same projection technique 
allows one to derive renormalization group equations. 

The projection technique is based on a clever insertion 
of projection operators into the equation of motion (|^). 
A projection operator is any operator V which satisfies 
= V; its complement, which is also a projection oper- 
ator, is denoted by Q := l — V. The projection operators, 
like the Liouvillian, are so-called "superoperators:" they 
do not act in Hilbert space, but in the space of observ- 
ables ("Liouville space"). For our purposes we consider 
projectors which project arbitrary vectors in Liouville 
space onto the subspace spanned by the unit operator 
and by the relevant observables {Ga}; i.e., for which 



'PA = A ^ Ae span{l, Ga} 



(5) 



For simplicity we assume that the Hamiltonian and hence 
the Liouvillian, as well as the relevant observables, are 
not explicitly time-dependent. In contrast, we allow the 
projector to depend on the expectation values {ga{t)} of 
the relevant observables, thus making it an implicit func- 
tion of time: V{t) = Vlgait)]; with the sole restriction 
that for any observable A, 



(pit)\^Vit)A)^0 



(6) 



For the time being we admit any projector which satisfies 
the two constraints (||) and (^). Only later, in Sec. iVB, 
we shall make a specific choice for Vit). 

Now let T{t', t) be the (super)operator defined by the 
differential equation 



d_ 

dt- 



-r{t' ,t) = -iQ{t')CQ{t')r{t' ,t) 



(7) 



with initial condition T{t,t) ~ 1. It may be pictured as 
describing the evolution of the system's irrelevant degrees 
of freedom. With its help the equation of motion for the 
selected expectation values {gait)} can be cast into the 
form 
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gait) = l{p{t)\V[t)CGa) 

~ f dt'{p{t')\V{t')CQ{t')T{t',t)Q{t)CGa) 
Jo 

+ iipmQiO)T{0,t)Q{t)CGa) , (8) 

for any time t > 0. Comparing this form of the equation 
of motion with the original form (^) we notice that, apart 
from the replacement (p| — > {p\V, there are two terms 
which are qualitatively new: (i) an integral ( "memory" ) 
term, containing contributions from all times between the 
initial and the present time; and (ii) a "residual force" 
term describing the effect of irrelevant components in the 
initial state. The physical meaning of both terms can be 
easily discerned if read from left to right: (i) At time 
t' < t relevant degrees of freedom (projected out by V) 
couple via an interaction (£) to irrelevant degrees of free- 
dom (projected out by Q), which subsequently evolve in 
time (T) and, due to a second interaction (C), acquire 
relevancy again, thus influencing the evolution of the rel- 
evant observable Ga at the present time t. (ii) Irrele- 
vant components in the initial state (Q) evolve in time 
(T) and, due to interaction (£), acquire relevancy at the 
present time t. 

In many practical applications the irrelevant compo- 
nent of the initial state can be shown to vanish, or else 
to be negligible; in which case the residual force term can 
be dropped from the equation of motion. What remains, 
is then the desired closed system of (possibly nonlinear) 
coupled integro-differential equations for the selected ex- 
pectation values {ga{t)}. The principal feature of these 
closed equations is that they are non-Mar kovian: future 
expectation values of the selected observables are pre- 
dicted not just on the basis of their present values, but 
based on their entire history. 

III. APPROXIMATIONS 

A. Second order perturbation theory 

Often the Liouvillian can be split into a free part and 
an interaction part, 

£ = £(0) + V , (9) 

corresponding to a decomposition H = iJ^'^^ -I- of the 
Hamiltonian. Provided free evolution does not mix rele- 
vant and irrelevant degrees of freedom, i.e., provided 

[c^"\r{t)]^o , (10) 

then in the memory term VCQ = WQ and QCP — 
QVV\ hence the memory term is at least of second order 
in the interaction. To second order, therefore, one can 
simply replace 

Q{t')T{t\t)Q{t) ^ Q{t')T^''\t\t)Q{t) , (11) 



where T*^") is a time-ordered exponential of QC^^^Q. Us- 
ing (^ and 2(^2)2(^1) = 2(^2), aU the 2's appearing in 
T(°) can be shuffled to the left and absorbed into Q{t')^ 
allowing one to replace further 

r(o) [f, t) ^ (i', t) := exp[i/:W • [t - t')] . (12) 

One thus obtains (without residual force) 

ga{t)^i{p{t)\V[t)CGa) 




dr {p{t - T)\V{t - T)VQ{t - t)U^'>^ (0, T)VGa) ■ 



(13) 



B. Markovian limit 

We have seen that predictions of future expectation 
values of the selected observables generally depend in a 
complicated manner on both their present expectation 
values and their past history. There are thus two distinct 
time scales: (i) the scale Trci -or several scales {r'^i}- 
on which the selected expectation values {ga{t)} evolve; 
and (ii) the memory time Tmom: which characterizes the 
length of the time interval that contributes significantly 
to the memory integral. Loosely speaking, the memory 
time determines how far back into the past one has to 
reach in order to make predictions for the further evo- 
lution of the selected observables. If this memory time 
is small compared to the typical time scale on which the 
selected observables evolve, r-a^cm ^ '''rci, then memory 
effects can be neglected and predictions for the selected 
observables be based solely on their present values. One 
may then assume that in the memory term gait') ~ ga(t) 
and hence replace 

V[ga{t')]^V[ga{t)] , 

ip{t')\r{t')^{p{t)\r{t) . (14) 

This is the Markovian limit. 

Closely related to the Markovian limit is the quasis- 
tationary limit: at times t ^ Tmam it no longer matters 
for the dynamics when exactly the evolution started, and 
hence in Eq. (j^ the integration over the system's history 
may just as well extend from —00 to t, rather than from 
to t. 

In the Markovian and quasistationary limits the equa- 
tion of motion -to second order and without residual 
force- simplifies to 

ga{t) ^ t{pit)\rit)CGa) 

- / dTip{t)\V{t)VQ{t)U^°'>{0,T)VGa). (15) 
Jq 

This approximate transport equation for the {ga(t)} shall 
be the basis of our further investigations. 
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IV. ANALYSIS 



p,,i{t):^Z{t)-^exp[-X-{t)Ga] , (23) 



A. Reformulation of the transport equation 

Starting from the approximate transport equation 
( |l5| ) , we eventually wish to discern two of the main fea- 
tures of macroscopic transport: (i) dissipation and (ii) 
the modification ( "renormalization" ) of the effective in- 
teraction. The latter will then, in Section lead on to 
the consideration of renormalization group equations. 

To this end we first split the free evolution operator 
into its symmetric and antisymmetric parts, 



with partition function 



+ i[Z^(°)(0,T)-Z^W(r,0)] 



use 



dr [U'-°^ (0, r) + U'-°^ (r, 0)] = 2T:d{C 



(o)^ 



(16) 



(17) 



in the symmetric part, and in the antisymmetric part 
exploit the liberty -thanks to the Markovian limit- to 
insert free evolution operators ZY''''(0,r) and ZY^°^(r, 0) 
in front oi'P{t) or Ga, respectively. We thus obtain 

gait) ^ iip{t)\Vit)Cesit)Ga) 

-7r{p{t)\V{t)VQ{mC^°^)VGa) , (18) 



where 



Cait) = C + SC{t) 



(19) 



denotes a (possibly time-dependent) "effective" Liouvil- 
lian determined by 



-[P{t)V{T)V{t),V{t)VVm , (20) 



and 



V{T)^h-'[V{T),*] (21) 

is the commutator with the interaction-picture operator 

Vir) ■.^U^"\0,t)V . (22) 

Next, we evaluate the reformulated transport equation 
( |l8| ) by making a suitable choice for the yet undetermined 
projector V{t). 

B. Robertson projector 

Associated with the expectation values (|]) of the se- 
lected observables is, at each time t, a generalized canon- 
ical state 



Z{t) :=tr cxp[-X''{t)Ga] 



(24) 



and the Lagrange parameters {A°(t)} adjusted such as 
to satisfy the constraints 



tr[prclit)Ga] = ga{t) 



(25) 



Here we have used Einstein's convention: repeated upper 
and lower indices are to be summed over. The general- 
ized canonical state p,-ci(t), among all states which satisfy 
the constraints (25), is the one which maximizes the von 
Neumann entropy 



S[p] := — fctr(plnp) 



(26) 



For this reason it may be considered "least biased" or 
"maximally non-committal" with regard to the unmoni- 
tored degrees of freedom; it is sometimes called the "rel- 
evant part" of the full statistical operator p{t). 

There exists a unique time-dependent projector Vnit) 
which projects arbitrary vectors in Liouville space onto 
the subspace spanned by the unit operator and by the rel- 
evant observables {Ga}, the projection being orthogonal 
with respect to the time-dependent scalar product 

{A;B)^'^ -.^ rdptY[p,,,{trA^p,,i{ty-^B] . (27) 



This projector Vv{,(t) satisfies both conditions (^) and (|^) 
and, moreover, can be shown to yield 



{pit)\Vnit) = (p,.el(i)| 



(28) 



at all times. This special choice for V{t), originally pro- 
posed by Robertsonj3'c2l has the important advantage 
that in contrast to other frequently used projectors such 
as the Mori projectoitl it permits the derivation of closed 
transport equations valid arbitrarily far from equilib- 
rium. We shall use this projector throughout the remain- 
der of the paper (and for brevity, we shall immediately 
drop the subscript 'R'). 

With the Robertson projector the transport equation 
( |l8| ) takes the form 

9a{t) = l{prelit)\C,sit)Ga) + Af,,(i)A^(t) , (29) 

where 

M,a{t) := TT{Qit)VGc; <5(£W)Q(i)VGa)(*^ (30) 

is a matrix whose eigenvalues are all real and non- 
negative. In this formulation it is particularly easy to 
distinguish the dissipative and non-dissipative parts of 
the dynamics: as we will show, the second term in ( p9| ) is 
solely responsible for dissipation; whereas the first term 
yields non-dissipative dynamics governed by a modified 
( "renormalized" ) effective Hamiltonian. 



4 



C. Dissipation 

An observer who monitors only the selected degrees of 
freedom does not have complete information about the 
system's microstate. A suitable measure for this lack of 
information is the entropy S[prei{t)] associated with the 
relevant part of the statistical operator. It is sometimes 
called the "relevant entropy." This relevant entropy gen- 
erally varies in time: it changes at a rate 



(31) 



Within our approximations -perturbation theory and 
Markovian limit- this rate may be evaluated by inserting 
the transport equation (p9|). Its first term does not con- 
tribute to the change of relevant entropy; only its second 
term yields a nontrivial contribution 



5[prelW] = fcAf,a(t)A=(t)A"(t) > 



(32) 



The relevant entropy thus increases monotonically, re- 
flecting dissipation and irreversibility of the macroscopic 
dynamics. It stays constant if and only if the second term 
in (29) vanishes ("adiabatic limit"). 

Of course, our finding represents one particular case 
of the more general i?-theorem. That the relevant en- 
tropy can never decrease, is a direct consequence of the 
Markovian limit and hence holds true whenever the sys- 
tem exhibits a clear separation of time scales.E3 



D. Effective Hamiltonian 

The non-dissipative part of the macroscopic dynamics 
is encoded entirely in the first term of (p9|). So in the 
adiabatic limit, which we shall consider from now on, the 
transport equation simplifies to 



ffa'^^(0 =»(Prcl(i)l'Ccff(i)Ga) 



(33) 



This is very similar to the original equation of motion 
(^, yet with p replaced by prei, and £ replaced by £eff- 
Here we make a special choice for the selected observ- 
ables, one that will directly lead to the renormalization 
group. We presume that we are interested in features of 
the macroscopic system (for example, its long-wavelength 
properties) which can be represented by observables act- 
ing merely in some subspace of the original Hilbert space 
(e. g., in the subspace spanned by all many-particle 
states with momenta below a given cutoff). Let the op- 
erator which projects the original Hilbert space onto this 
selected subspace be denoted by P, and its complement 
by Q = 1 — P. Selected observables are then all those 
of the form PAP, with A being an arbitrary Hermitian 
operator. This choice of relevant observables gives rise 
to a particularly simple representation of the Robertson 
projector. 



VA = PAP 



tviQA) 
trQ 



Q 



(34) 



which no longer varies in time. 

We now decompose the microscopic Hamiltonian H 
into a non-mixing "free" part 



(35) 



i/(0) = pjjp QJJQ 

and a mixing "interaction" 

V = PHQ + QHP , (36) 

and split the Liouvillian correspondingly. Since 

[C^°\V]^0 , (37) 

this decomposition is suitable for perturbation theory. 
With the property 

VVV = (38) 
we can immediately evaluate Eq. ( pO[ ) to obtain 

d£ = n-^[dH,*] , (39) 
where 6H is given by 



SH ■ 



2h 



dr [V, V{t)] 



(40) 



This yields then 

V£esVA^n~^[H,s,PAP] VA , (41) 

where we have identified the effective Hamiltonian 

iJcff = PHP + E (42) 

which is not just the projection PHP of the original 
Hamiltonian, but contains an extra term 

^--^{PHQ^QHP + h.c. 

+ ^{iPHQS{C^°^) QHP + h.c.} (43) 

stemming from SH. In S the contributions which involve 
(5(£'-*'-' ) may be omitted as long as the P- and Q-scctors of 
Hilbert space are associated with clearly distinct energies. 

Our result for the effective Hamiltonian, which we have 
obtained within the general framework of transport the- 
ory, is very similar to the Block- Feshbach formula known 
in the theory of ruiclear dynamics,Eil or to Anderson's 
poor man's scaling.i3 Below we wish to demonstrate how 
this result can be utilized to derive renormalization group 
equations for a variety of physical systems. 
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V. RENORMALIZATION GROUP FOR 
INTERACTING QUANTUM GASES 

A. Hamiltonian and Ground State 

As an illustration of the above general result we shall 
investigate the low-temperature properties of interacting 
quantum gases, i.e., the effective dynamics of low-energy 
excitations above their many-particle ground state. We 
assume the microscopic dynamics of the gas to be gov- 
erned by a Hamiltonian of the form 

^^^k ■.alak:+'^^{lk\V\ji)± :a|4ajaj: , (44) 

k ijkl 

with kinetic energy i?kin and a two-body interaction 
V2^2- The single-particle energies e/c include the chemi- 
cal potential. Annihilation and creation operators obey 
[ai,aj] ip — Sij for bosons (upper sign) or fermions (lower 
sign), respectively. 

Each term in the Hamiltonian is normal ordered 
(: ... :) with respect to the nonintcracting many-particle 
ground state. For bosons this ground state has all parti- 
cles in the lowest energy, zero momentum single-particle 
mode; while for fermions it consists of a filled Fermi sea 
with all momentum modes occupied up to some Fermi 
momentum Kp. (For simplicity, the Fermi surface will 
be taken to be spherical.) The explicit normal ordering 
of the Hamiltonian is redundant in the bosonic case. In 
the fcrmionic case, on the other hand, it means shuffling 
all operators which annihilate the fermionic vacuum (a^ 
for states above the Fermi surface, a| for states below 
the Fermi surface) to the right, all others (a|, for states 
above, Uk for states below the Fermi surface) to the left, 
thereby changing sign depending on the degree of the 
permutation. 

We will assume that, at least to a good approximation, 
the essential features of the ground state survive even in 
the presence of interaction. More specifically, we will as- 
sume that in the case of interacting bosons the ground 
state still has most particles in modes with zero, or at 
least very small, momentum; and that in the case of in- 
teracting fermions there still exists a well-defined Fermi 
surface. Low-energy excitations then correspond to the 
promotion of bosons from small to some slightly higher 
momentum, or of fenaions from just below the Fermi sur- 
face to just above it.EHl At low temperature the regions 
of interest in momentum space are therefore the vicinity 
of the origin (bosons) or the vicinity of the Fermi surface 
(fermions) , respectively. 

B. Mode Elimination 

We now wish to devise a systematic procedure for fo- 
cusing onto these regions of interest. To this end we con- 



sider effective theories (i) in the bosonic case for modes 
within a sphere around the origin, of radius A; and (ii) 
in the fermionic case for modes within a shell inclosing 
the Fermi surface, of mean radius Kp and thickness 2A 
(where A <C Kp). Whereas in the limit of large A one re- 
covers the original, full theory, the opposite limit A ^ 
yields the desired low-energy effective theory. In order 
to interpolate between these two limits we proceed in in- 
finitesimal steps: we lower the cutoff from some given 
A(s) to 

A(s -I- As) := exp(-As)A(s) , As > (45) 

with As infinitesimal, thereby discarding from the the- 
ory momentum modes pertaining to an infinitesimal shell 
(in the fermionic case: two shells) of thickness AA = 
A(s) As; then determine the effective dynamics of the 
remaining modes; then eliminate the next infinitesimal 
shell, again determine the effective dynamics of the re- 
maining modes, and so on. After each infinitesimal step 
we obtain a new effective Hamiltonian with slightly mod- 
ified coupling constants. These may also include novel 
couplings which had not been present in the original the- 
ory — in fact, the mode elimination procedure will typi- 
cally generate an infinite number of such novel couplings. 
But in many cases only a few coupling constants will 
change appreciably and thus suffice to study the physi- 
cal system at hand. How these coupling constants evolve 
as the flow parameter s increases and hence the cutoff 
A(s) approaches zero, can then be described by a small 
set of coupled differential equations. Modulo trivial scal- 
ing, these are the renormalization group equations of the 
theory. 

At a given cutoff A the many-particle Hilbert space 
(Fock space) for bosons is spanned by the particle-free 
vacuum |0b) and all n-particle states (n = 1 . . . oo) 

n 

|ki...k„) cx[]at(kO|Ob) , |k,|<A , (46) 

i=l 

where the {k^} denote the particle momenta and {a^(ki)} 
the associated bosonic creation operators. 

The fermion Fock space, on the other hand, is spanned 
by the filled Fermi sea (fermionic vacuum) |0f) and all its 
excitations which have particles above the Fermi surface 
and/or missing particles ("holes") below it, all within 
a shell of thickness 2A. In order to cast this into a 
mathematical formulation it is convenient to change co- 
ordinates, from the true particle momenta {K^} to little 
( "quasiparticle" ) momenta 

k, := (|K,|-i^;^)K, (47) 

and additional discrete labels 

a, :=sign(|K,|-ifj.) . (48) 

This coordinate transformation K (k, a) is invertible 
except for modes which lie exactly on the Fermi surface. 
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States above the Fermi surface are labeled cr = 1, while 
those below are labeled a = —1. The allowed excitations 
in fermion Fock space then have the form {n — 1 . . .00) 



|k± . . . k±) « Y[[eia,)a^ik,,a,) + 0(-a,Oa(-k„ a,)]|Of) , 

|k.|<A , (49) 



where the {k^} denote the momenta of particles (+) 
or holes (— ), respectively, and {a^, {a} the associated 
fermionic creation and annihilation operators. For sim- 
plicity, we have omitted any spin quantum numbers. 

It is now obvious which form the projection operator 
will have that is associated with the infinitesimal cutoff 
reduction (^5|): if applied to any of the excitations ( ^ ) 
or (^9|) it will simply multiply the respective state by a 

product of 9 functions, Jli ^(^ ~ '^'^^m-''^^ I): to enforce 
the new cutoff constraint. 



C. Modification of the Two-Body Interaction 

Each mode elimination will yield an effective Hamilto- 
nian that will generally contain a slightly altered mass, 
chemical potential, two-body interaction, etc., and pos- 
sibly new interactions such as an effective three-body in- 
teraction. Here we shall restrict our attention to the 
modification of the two-body interaction. This modifica- 
tion is entirely due to the extra term E (Eq. (^)) in the 
effective Hamiltonian, 

^ = I E Y.(^k\Vm±{dc\V\ba)± 

I abed ijkl 

X P -.alalajai-. Q-^^Q -.al^alabaa- P + h.c.^ , (50) 

where to the given order in perturbation theory C^^'> just 
coincides with £kin- The two projectors P at both ends 
of the operator product ensure that all external momenta 
lie below the new, reduced cutoff; whereas the projectors 
Q in the center force at least one internal momentum to 
lie in that infinitesimal shell which has just been elim- 
inated. Therefore, at least one pair of field operators 
must pertain to the eliminated Q-modes, and hence be 
contracted. The product of the remaining six field op- 
erators can then be rearranged with the help of Wick's 
theorem to yield a decomposition 



(51) 



each E„ being a normal ordered product of n field op- 
erators. The various terms shift the ground state en- 
ergy (n = 0), modify the mass, the chemical potential, 
or more generally the form of the single-particle disper- 
sion relation (n = 2), modify the two-body interaction 



(n = 4), and generate a new effective three-body inter- 
action (n = 6). 

As we want to focus on the modification of the two- 
body interaction, we consider only the term with n — -i. 
Neglecting the energy of the external modes we find for 
bosons {h = 1) 



A(lk\V\ji)+ = -A 



and for fermions 



E or. L A mV\ba) + {ab\V\jz) + 

.ah 

(52) 



A(;fc|T/|ji)_ = A 



with 



^(zs) _ 

lk\ji 



ah 



(ZS) ^(ZS') ^(BCS) 

lk\ji i^lji If^lj'i 



(53) 



its cross term 



x{la\V\hi)-{hk\V\ja). 



^(ZS') ^ _^(ZS) 
lk\ji 



(54) 



(55) 



and 



^(BCS) _ _^rV- 

^ik\ji - ^\-2^ 



ah 



2(ea + Eb) 



x{lk\V\ba)^{ab\V\ji)-] 



(56) 



The A in front of the sums signifies that at least one 
of the internal modes (a, b) must lie in the eliminated 
shell. In the bosonic case the modification of the two- 
body interaction can be associated with a 1-loop "ladder" 
diagram. In the fermionic case, on the other hand, there 
are three distinct contributions which, with hindsight, 
may be identified with "zero sound" (ZS,ZS') and BCS 
diagramsJlil The ZS contribution and its cross term ZS' 
account for particle- hole excitations (cto — ±1, cTf, = Tl), 
while the BCS term describes 2-particle (tTa = crb = +1) 
or 2-hole (da = ai, — —1) excitations. 

The above formulae serve as the starting point for the 
investigation of a variety of specific physical systems. 
With their help one can derive such diverse results as 
the 1-loop renormalization group equation for an inter- 
acting Bose gas, the 1-loop /3 function of 0"* theory, the 
screening of fermion-fermion interactions, or the BCS in- 
stability. Details of these applications are presented in 
the following section. 



VI. EXAMPLES 

A. Bosons witli point interaction 

For spinless bosons with point interaction (5 function 
potential in real space) it is 
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2U 

{lk\V\ji)+ = -^(5k,+k,.k,+k, 



(57) 



with the Kronecker symbol enforcing momentum con- 
servation, being the spatial volume, and U the cou- 
pling constant. Provided the magnitude of the external 
momenta k^, k^ is negligible compared to the cutoff A, 
momentum conservation implies that the internal modes 
a, b must both lie in the eliminated shell, and that hence 
— ~ ^K- Application of the general formula ( |5^ ) 
then yields 



where the sum 



|k„Mk5|<A 



k„+k(, 



(58) 



A 



(5k;+kj,k„+ki, 

|k„|,|k5|<A 



A '5k„-k, 
|ka|.|k6|<A 

E 1 

|k„|(E[A-AA,A] 



(59) 



simply counts the number of eliminated states. For a 
spherical cut in momentum space this number of states 
is given by 



E 

|k„|e[A-AA,A] 



1 



(60) 



with p(eA) denoting the density of states at the cutoff. 
With AA = A As we thus obtain the flow equation 



AC/ = 



dlnA 2n 



(61) 



For a dilute gas of nonrelativistic bosons in three spa- 
tial dimensions, with mass m, dispersion relation ea = 
A^/2m and density of states p{eA) = ^lmA/2Tr^ the flow 
equation reduces to 



AC/ = 



mA 
2^ 



U'-As 



(62) 



By its very definition the sequence of effective theories 
retains complete information about the system's low- 
energy dynamics. Observables pertaining to this low- 
energy dynamics are therefore unaffected by the succes- 
sive mode elimination, and hence independent of s. For 
example, the scattering lengthen 



m 
47r 



U{s) - U{s)' 



|p|<A(s) 



cPp m 
(27r)3 ^ 



stays constant under the flow (|62|), as the s-dependence 
of the parameters U and A just cancels out (up to third 
order corrections). 



B. The link to 0* theory 

There is an interesting relationship between the result 
( pi] ) and the 1-loop /3 function for real 0* theory. The 
(]r Hamiltonian describes the dynamics of coupled anhar- 
monic oscillators. It reads, in three spatial dimensions, 
H = _ff + V with kinetic energy 



k«k 



(64) 



and interaction 



in E n7i^(«'^^+«^-kJ^Sk.o- (65) 

4 a=l ^ 



kik2k3k 



Here m denotes the mass, the spatial volume, g the 
coupling constant, and ek the single-particle energy 



ek 



(66) 



The field (j) its conjugate momentum tt are time- 
independent (Schrodinger picture) operators which sat- 
isfy the commutation relations for bosons, and a, 
are the associated annihilation and creation operators. 
While the kinetic part of the Hamiltonian is normal or- 
dered (: ... :), the interaction is not. 

When expressed in terms of annihilation and creation 
operators the Hamiltonian takes on a form which is very 
similar to that of the quantum gas Hamiltonian (|44|). 
More precisely, the (f)^ Hamiltonian contains a Bose gas 
Hamiltonian with two-body interaction matrix element 



{ik\vm 



4^ g 

2) iin^ye-eje^' 



'k,+k,-,kfc+k. 



(67) 



The derivation of a flow equation for g can now proceed 
in the same vein as that for U, again starting from Eq. 
(^^ . Now, however, apart from 2 — > 2 particle scattering 
the (p'^ Hamiltonian with its novel interactions a^a^a^a, 
a^a^a''a'' etc. also permits 2 — > 4 and 2^6 scattering. 
Therefore, in Eq. ( p2[ ) the intermediate state may be not 
just \ab), but also \abik), \abil), \abjk), \abjl) or \abijkl). 
As long as the magnitude of the external momenta is 
negligible compared to the cutoff, it is in all six cases 
fa = £6 = EA and 



(63) {lk\V\ 



\vm+ = 



{lk\V\jt)+S^ 



(68) 



Hence in order to account for the larger set of allowed 
intermediate states we merely have to introduce an extra 
factor 6, and obtain thus 
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dln £A 3p(£a) ^2 
din A Srtel 



(69) 



For A m it is CA = A, p{eA) — fle\/2TT^, and the flow 
equation reduces to 



(70) 



in agreement with th e wel l-known 1-loop result for the /3 
function of t/)^ theory.cJO 



C. Screening of fermion-fermion interactions 

We consider nonrelativistic fermions in spatial dimen- 
sion d {d > 2) which interact through a two-body inter- 
action 



x6 



(71) 



duly antisymmetrized to account for Fermi statistics, and 
with {sq} denoting the spin quantum numbers and q, q' 
the respective momentum transfers 



q 

q': 



Kfe 



Kfe 
K, 



We investigate scattering processes for which 
0< |q|,A« |q'|,|K, + K,|,if^^ . 



(72) 



(73) 



In this regime only the ZS contribution (|5J) can signif- 
icantly modify the two-body interaction; its cross term 
ZS' (Eq. (H)), as well as the BCS contribution (|56|), 
are suppressed by a factor A/Kp. This can be seen 
directly from the geometry of the Fermi surface. The 
three constraints on the intermediate state: (i) both 
Ka and Kf, lie in the cutoff shell of thickness 2A; (ii) 
more stringently, one of them lie in the infinitesimal shell 
to be eliminated; and (iii) — K{, = — q' (ZS') or 
Ka + Kh = Ki + Kj (BCS), respectively — reduce the 
momentum space volume available to the internal mo- 
mentum Ka to 0(A'^~2AAA). In contrast, for |q| ^ A 
the ZS contribution with its condition Ka — Kb = — q 
allows a momentum space volume of the order A'^^^AA. 

To evaluate the ZS contribution at some given momen- 
tum transfer q, we first define the angle i) between — q 
and the internal momentum Kq, 



cos ■& 



Ka 



|q||Ka| 



(74) 



change coordinates from original (K) to little (k) mo- 
menta, and write, up to corrections of order \q\/Kp, 

ia- (-b^ fFdKol - |Kb|) = VF{<Ja\'ka\ - (Jb\'i^b\) 

^VF\q\z (75) 



with vp denoting the Fermi velocity. Next we note 
that the term with O{aa)0{—(Jb) and the term with 
d{—(Ta)d{ab) yield identical contributions; therefore it 
suffices to consider only the first term and then multi- 
ply it by two. Finally, assuming that in the interaction 
matrix element (|7l|) the cross term is negligible, 



inq')l«inq)l 



(76) 



the two matrix elements in Eq. ( p4]) can simply be re- 
placed by V{q)^ modulo Kronecker symbols for spin and 
momentum conservation. By virtue of these Kronecker 
symbols one of the two summations over internal modes 
collapses trivially, leaving 



AV^(q) = A 

VF 



5]0(A~|ka|)0(|ka|-|q|z-KA) 
,0(aa)0(|q|z-|ka|)" 



V{ciY . (77) 



Here the first two 9 functions explicitly enforce the sharp 
cutoff constraint for both ka and kf, (|ka|, |kf,| < A), while 
the latter two functions enforce Ca = 1 and CTf, = —1, 
respectively. Under these constraints it is always ■& G 
[0,7r/2) and hence z > 0. 

The above equation can be immediately integrated 
from cutoff A ^ |q|, |ka| (symbolically, A — > oo) down 
to A <C |q|, jkal (symbolically, A ^ 0), to yield the total 
modification of the two-body interaction: 



1 



1 



Kfr(q) Vbare(q) 
VF ^ 



|q|z-f A)0(aa)^(|q|z~|ka|) 



|q|^ 





(78) 



At the lower bound (A 0) the various conditions im- 
posed by the Q functions cannot all be satisfied simultane- 
ously, and therefore the product of Q functions vanishes. 
At the upper bound (A — > oo), on the other hand, the 
cutoff constraints imposed by the first two 9 functions 
are trivially satisfied and thus can be omitted. In this 
case the sum over a is evaluated by turning it into two 
integrals, one over a radial variable such as |ka| or ea, 
the other over the solid angle. At a given solid angle and 
hence given z, the fourth function restricts the radial 
integration to the range |ka| € [0, |q|z] or, equivalently, 
ta S [0; ^^-F|q|^]- This energy interval in turn corresponds 
to a number [p(eir)wi?|q|z] of states, p(tF) being the den- 
sity of states at the Fermi surface. (It includes a factor to 
account for the spin degeneracy.) The integration over 
the solid angle is constrained to a semi-sphere, due to 
•d e [0, 7r/2), and hence reduced by a factor 1/2 as com- 
pared to a full-sphere integration. Altogether we obtain 



1 



1 



(79) 
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and thus 



VII. CONCLUSION 



Kff(q) = 



1 



Vbarc(q; 



+ p{f-p) 



(80) 



This result describes iie weh- known screening of fermion- 
fermion interactions .c3 



D. BCS instability 

Our last example pertains to fermions with an attrac- 
tive pairing interaction 



(81) 

which is the simplest form of BCS theory.0 Due to the 
pairing condition Kj — — K^, K; = — Kfc it is impossible 
to satisfy in the ZS and ZS' terms the requirement that at 
least one of the internal modes be in the eliminated shell. 
Hence only the BCS term ( |56| ) can modify the coupling 
constant. In the BCS term there are contributions with 
9{(Ta)S{<7b) and 9{—(Ja)0{—ah), respectively, which yield 
identical results; therefore it suffices to consider only the 
first contribution and then multiply it by two. The pair- 
ing condition implies ea = £6 = ^a, which for modes in the 
upper {(J = +1) eliminated shell is given by ea = wfA. 
The eliminated shell itself covers an infinitesimal energy 
interval of width [uj^AAs], which in turn corresponds to 
a number [p(ei?)ui?AAs] of states. Of the two summa- 
tions over internal modes one collapses trivially due to 
momentum and spin conservation, leaving 



AV 



2WFA 



-A 



(82) 



From this flow equation for the BCS coupling V{s) we 
immediately conclude that as long as the initial coupling 
V{Q) is positive, V{s) diverges as s ^ oo. This indicates 
the occurence of binding ("Cooper pairs") at very low 
temperatures. Furthermore, we can again convince our- 
selves that the sequence of effective theories retains com- 
plete information about the system's low-energy dynam- 
ics: Jew-energy observables such as the zero-temperature 
gapcZl 



Aq = 2A(s) exp 



p{^f)V{s) 



(83) 



do not depend on the flow parameter s and are thus un- 
affected by the successive mode elimination. 



We have succeeded in linking renormalization to trans- 
port theory. Our line of argument proceeded from the 
exact microscopic dynamics, via the Nakajima-Zwanzig 
projection technique, to a macroscopic transport equa- 
tion for selected expectation values, and then, via sec- 
ond order perturbation theory, Markovian limit, choice 
of the Robertson projector and a suitable rearrangement 
of terms in the transport equation, to the approximate 
effective Hamiltonian which governs the non-dissipative 
part of the macroscopic dynamics. We investigated the 
ramifications of this result for the low-energy dynamics of 
interacting quantum gases: contracting the set of selected 
expectation values by discarding iteratively those observ- 
ables which pertain to short- wavelength excitations, we 
obtained a sequence of effective Hamiltonians which de- 
scribe the dynamics on successively larger length scales. 
We then focused on the two-body interaction in these 
effective Hamiltonians and convinced ourselves -in sev- 
eral rather diverse applications- that it varies in accor- 
dance with 1-loop renormalization group equations. We 
have thus substantiated our original claim that renor- 
malization group equations can be obtained within the 
Nakajima-Zwanzig projection approach, and that hence 
renormalization can be embedded into the general math- 
ematical framework of transport theory. 

There remain many open questions worth investigat- 
ing. Clearly, a unified theoretical framework for dissi- 
pative transport and renormalization will have to prove 
its merits in new applications where the conventional ap- 
proaches fail: for example, when time scales are no longer 
well separated and the Markovian limit ceases to be justi- 
fied; or when in the course of successive mode elimination 
one starts to discard states with a finite population, thus 
introducing dissipation into the effective macroscopic dy- 
namics. It is my hope that the present paper will help 
stimulate research efforts in these and related directions. 
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